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On the pulse-width statistics in radio pulsars 

Marcin Kolonko^'^, Janusz Gil^ and Krzysztof Maciesiak^ 

ABSTRACT 



The Monte Carlo simulations of pulsar periods, pulse-widths and magnetic 
inclination angles are performed. Using the available observational data sets we 
study a possible trial parent distribution functions by means of the Kolmogorov- 
Smirnov significance tests. We also use an additional condition that the numbers 
of generated interpulses, whether from both magnetic poles or from single pole, 
are at the observed levels. We conclude that the parent distribution function of 
Q>^ ' magnetic inclination angles is neither flat nor cosine but it is a more complicated 

■ function with a local maximum near ol = 25° and another weaker one near 

CS| ■ a = 90°. The plausible distribution function of pulsar periods is represented by 



the gamma function. The beaming fraction describing the fraction of observable 



O ■ radio pulsars is about 0.12. 



O , 

^ , 1. Introduction 



Statistical studies of the pulse-width in mean profiles of radio pulsars are an important 



H ' tool for investigations of the pulsar radiation geometry. One especially important parameter 
that can be derived from such studies is the inclination angle between the magnetic and 
the spin pulsar axes. Early studies were carried out by Henry & Paik (1969); Roberts & 
Sturrock (1972, 1973); Backer (1976) and Manchester & Lyne (1977). Since the amount of 
the available data was small, these papers suffered from problems of small number statistics. 
A more complete work was performed by Proszyhski (1979) and Lyne & Manchester (1988), 
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who analyzed samples of about 200 pulse-width data measured near 400 MHz. Although the 
database used in these papers was quite rich, the pulse-width measurements were contami- 
nated by the interstellar scattering dominating at low radio frequencies. More recently Gil & 
Han (1996, GH96 hereafter) compiled a new database of 242 pulse-widths Ww (correspond- 
ing to about 10% of the maximum intensity) measured at a higher radio frequency (near 1.4 
GHz) , which was relatively unbiased as compared with the lower frequency data. GH96 used 
their pulse-width database to perform Monte Carlo simulations in an attempt to derive the 
distribution statistics of pulsar periods, pulse widths, magnetic inclination angles and rates 
of the interpulse occurrence. They concluded by comparing the simulated and observed (or 
observationally derived) quantities that the observed distribution of the inclination angles 
resembles the sine function following from the flat (random) distribution in the parent pop- 
ulation, and that the probability (beaming fraction) of observing a pulsar was about 0.16. 
GH96 also pointed out that the rates of interpulse occurrence should be considered as an 
important aspect of pulsar population studies. 

On the other hand, Tauris & Manchester (1998, henceforth TM98) using a differ- 
ent method based on an analysis of the indirectly derived polarization position angles and 
magnetic inclination angles concluded that the observed distribution of the latter is cosine- 
hke rather than the sine-like suggested by GH96. They also obtained the beaming factor 
0.1 ±0.02, considerably lower than 0.16 obtained by GH96. TM98 pointed out a likely source 
of this discrepancy, namely the incorrect assumption used by GH96 that the observed dis- 
tribution and the parent distribution of pulsar periods are similar. Recently, Zhang, Jang, 
& Mei (2003, henceforth ZJM03) followed the Monte-Carlo simulation scheme developed by 
GH96. ZJM03 argued that both the parent distribution function and the observed distri- 
bution of pulsar periods can be modelled by the gamma function but with different values 
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of free parameters, and their Monte-Carlo simulations included searching for a 2-D grid of 
these parameters. As a result, ZJM03 concluded that indeed the cosine-like distribution 
(suggested by TM98) is a much more suitable to model the inclination angles in the parent 
pulsar population than the fiat distribution (suggested by GH96). They argued that a most 
plausible parent distribution is a modified cosine function, which has a peak around 25° and 
another weaker peak near 90°. They also obtained the beaming factor ~ 0.12, consistent 
with the result of TM98. 

As emphasized by ZJM03 in the conclusions of their paper, neither these authors nor 
TM98 considered potentially important constraints related to the observable interpulse emis- 
sion. Although GH96 did consider the issue of interpulse emission, but their estimate of rates 
of occurrence was not quite correct (see §3.4 in this paper). Moreover, their statistical anal- 
ysis was biased by incorrect assumption mentioned above concerning parent distribution of 
pulsar period. In this work we follow the simulation scheme of ZJM03 but we include the 
actual rates of interpulse occurrences. We demonstrate that a pure cosine distribution of 
the magnetic inclination angles generates much too few interpulses as compared with ob- 
servations and should be rejected as a plausible distribution function. We found that the 
modified cosine function of ZJM03 (see eq. [9]) is much better in this respect, mainly be- 
cause of the weak second maximum near a — 90°. Both the above distribution functions can 
reproduce the observed distributions of pulse-widths, pulsar periods and inclination angles 
almost equally well, and thus the observed interpulse statistics provides the most restrictive 
constraint discriminating between different trial distribution functions. 
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2. 



Basic formulae 



2.1. 



Pulse- widths 



It is generally accepted that the pulsar radio emission is relativistically beamed along 
the open dipolar field lines (to within I/7, where 7 ~ 100 is the Lorentz factor of the emitting 
sources). Thus, the pulse width Wio at the level of 10% of the maximum profile intensity 
can be written as 



(Gil 1981), where pio is the beam- width corresponding to 10% intensity level, a is the 
inclination angle of the magnetic axis to the spin axis, /3 is the impact angle of the closest 
approach of the hne-of-sight to the magnetic axis, and thus ^ = a + (3 is the observer angle 
between the spin axis and the line-of-sight. It is worth noting that equation (1) assumes 
the symmetry of a pulsar beam with respect to the fiducial plane containing both the spin 
and the magnetic pulsar axes. Thus, the appropriate pulse widths database should contain 
only pulsars with symmetrical profiles for which W — 2(fi, where (f is the pulse longitude 
measured from the fiducial phase corresponding to the fiducial plane (e.g. GH96). 



The opening angle of the radio beam (beam-width) corresponding to 10% intensity 
level can be derived from pulse width measurements Wiq (and a, j3 values derived from the 
polarization data) in the form of the so-called p — P relation. Lyne & Manchester (1988) 
obtained pio ~ 6.°5P~^/^ for the pulse-width data at 408 MHz, which scaled to 1.4 GHz 
writes 




(1) 



2.2. Opening angles 



Pio ~ 



5.°8P 



-0.33 



(2) 
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However, Biggs (1990) reanalyzed the same data sample and concluded that 



Pio ~ 



5.°6P 



i-0.5 



(3) 



Rankin (1993a,b) analyzed a large number of available pulse- width data interpolated 
to the frequency ~ 1 GHz and obtained a bimodal distribution of the opening angles. This 
result was later confirmed at frequency ~ 1.4 GHz by Gil, Kijak & Seiradakis (1993, GKS93 
henceforth) and independently by Kramer et al. (1994) and the resulting p — P relation reads 



with smaller angles preferred at shorter periods. GKS93 argued that 4.°9 is preferred in 
80% of pulsars with P < 0.7 s, and this constraint was used by GH96 (their Table 2), while 
ZJM03 used 4.°9 for P < 0.7 s and 6.°3 for other periods (their Table 1). We denote the 
former option by equation (4a) and the latter one by equation (4b) in Table 2. It is worth 
noticing that p — P relations expressed by equations (2)- (4) are, to some extent, equivalent. 
In fact, equation (2) seems to reflect the tendency of smaller angles being preferred at short 
periods in equation (4), and equation (3) represents the average of both values appearing 
in equation (4). Nevertheless, we use all the above relations in our simulation procedure 
(section 4). 



It is reasonable to assume that both the rotation axis and the observer's direction are 
randomly oriented in space. Thus, the probability density function for the observer angle 

^ = a + P is 




(4) 



2.3. 



Pulsar axes 



/(0=sin(a + /3). 



(5) 
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However, the distribution of the inchnation angle may depend on many unknown factors, and 
the most common probabihty density functions in the parent population of pulsars include 
the flat function 

/(«) = - (6) 

TT 

for the random distribution, the sine function 

f{a) — sina (7) 
for the sine-like distribution, and the cosine function 

f{a) — cos a (8) 

for the cosine-like distribution. One can consider some more complicated probability den- 
sity functions. For example, ZJM03 argued that the plausible parent distribution of the 
inclination angle can be described in the form 

A R 

f(a) = - \ =- (9) 

■'^ ^ cosh(3.5(a - 0.43)) cosh(4.0(a - 1.6)) ' 

where A — 0.6 and B — 0.15. This function has a weak local maximum near a — 7r/2, which 
is an important feature with regard to the interpulse occurrence statistics. We test all the 
above functions in our Monte Carlo simulation procedure (section 4). 

2.4. Pulsar periods 

GH96 demonstrated that the observed distribution of periods of 516 pulsars with 4.2 s > 
P > 0.05 s (Pulsar Catalog Taylor, Manchester & Lyne 1993) can be well fitted by the gamma 
function 

/(P) = Gox'^-'e-^ (10) 
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with X — P/m, where m — 0.3 and a — 2.5 (Go is the normahzation constant dependent 
on a. ZJM03 used this function to fit the period distribution in a much larger sample of 
1164 pulsar periods 4.2 s > P > 0.05 s and obtained m = 0.278 and a = 2.277. Although 
the parent distribution of periods is certainly different from the observed distribution, it is 
convenient to represent the former also in the form of a general gamma function (eq. [10]), 
with values of m and a treated as free model parameters that can be derived in 2-D grid 
search implemented into the Monte Carlo simulations (section 4). 

For a comparison one can also try some other trial probability density distribution 
functions, like the lorentzian distribution function 



l + (P-Xo)Va^' 



f{P) = TTT7r^h^2, (11) 



or even the gaussian distribution function 



V^TTCTo 



(12) 



2ao 

Again, the values of free parameters can be derived in 2-D grid search within the Monte 
Carlo simulations (section 4). 



2.5. Detection conditions 

Following the arguments given by Lyne & Manchester (1988); Gil & Han (1996) and 
Mitra & Deshpande (1999) we first assume that the pulsar beam is circular or almost circular. 
If this is the case, then the detection condition (see GH96 for details) is: 

P>\P\ (13) 

for the typical Main Pulse (MP) emission. However, one should also include a rare interpulse 
(IP) emission, seen in a few percent of pulsars. Two models were proposed to explain this 
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emission occurring about 180 degrees of longitude from the centroid of the MP: the double 
pole (DP-IP) model, in which MP and IP originate from the opposite magnetic poles (Rickett 
& Lyne 1968; Manchester & Lyne 1977), and the single pole (SP-IP), in which both MP and 
IP are associated with the same pole. The detection condition for DP-IP model is 



(e.g. GH96), which favors a nearly orthogonal rotators (a ~ 7r/2). Both MP and IP in 
pulsars that belong to DP-IP class show relatively narrow duty cycles, contrary to the broad 
profiles in the SP-IP class associated with nearly aligned rotators (a ~ 0). There are two 
possibihties within the SP-IP: in the first one the MP and the IP represent two cuts through 
one conical beam (Manchester 1977; Manchester & Lyne 1977). The other possibility is 
when the line-of-sight stays in the overall pulsar beam for the entire pulsar period, so the 
MP and IP correspond to cuts through two nested conical beams (Gil 1983). The detection 
conditions for these two versions of SP-IP models are 



respectively. Both these conditions are checked alternatively. One should mention that only 
the latter condition was used by GH96, which resulted in too low rates of SP-IP as compared 
to observations. 

One can also consider a quite natural possibility that the pulsar beam has a tendency 
to a meridional compression, with the ratio of minor to major ellipse axes 



p > TT — 2a — (3 



(14) 



(15) 



and 



p>2a + (3, 



(16) 




(17) 
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depending on the inclination angle (Biggs 1990; McKinnon 1993). Then the appropriate 

detection conditions (generalizations of equations [13]-[15]) are: 

< Vi?V + /5'(l-^')- (18) 

for the Main Pulse emission, 

rR> t: -2a- 13 (19) 

for the DP-IP emission, and 

rR>2a + (3 (20) 
or rR > \/2oP- + for the SP-IP emission, where 



^ = + (21) 

is the latitudinal beam dimension. Note that R< \ and r > where p is the observationally 
deduced beam radius (see GH96 for details). 

3. Observational data 

In our statistical analysis we compare the simulated distributions with the directly 
{WiQ, P, IP occurrences) or indirectly (a) observed data by means of the Kolmogorov- 
Smirnov (K-S henceforth) significance tests. Using the numerical methods given by Press et 
al. (1992) we compute the maximum distance T> (maximum difference between cumulative 
distribution functions corresponding to the observed and simulated data sets) as well as the 
significance V of any non zero value of V. The value of V represents the probability that 
both the simulated and the observed data sets are drawn from the same parent distribution. 
We adopt an arbitrary criterion that V > 0.1% for distributions of P, W^q and a at the 
same time. 
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3.1. Puls£ir periods 

Using the available database of pulsar parameters^ we selected a sample of 1165 periods 
P with 0.02 s < P < 8.52 s (Fig. 1). We rejected all recycled and binary pulsars, since they 
represent different period and magnetic inclination angle populations than typical pulsars 
(the recycled pulsar B1933+16 with P = 0.089 s was also excluded). This sample of the 
observed pulsar periods can be fitted by the gamma function expressed by equation (10) 
with m — 0.28 and a — 2.28 (see upper panel in Fig. 1). 

3.2. Pulse-widths 

We use the database of 238 pulse-width measurements Wio taken at the frequency of 
about 1.4 GHz (upper panel in Fig. 2), compiled by GH96 (their Table 1 and Fig. 1)^. These 
measurements were carefully selected from the available databases to satisfy all criterions 
imposed by the symmetry of equation (1), (see comments below this equation and GH96 for 
more details). ZMJ03 demonstrated that another available database (containing 265 pulsars) 
of Gould & Lyne (1998) is equivalent to that of GH96 in the sense that the values of Wio 
are roughly the same in both these databases. 

3.3. Inclination angles 

While the values of Wio and P as well as the interpulse occurrences are direct observa- 
tional quantities, the values of inclination angles a can only be indirectly observed by means 
of the polarization measurements (Lyne & Manchester 1988; Rankin 1990, 1993a,b; Gould 

^http: / /www.atut.csiro. all/research/ catalogue/ 

^We excluded 4 recycled and binary pulsars from the database of GH96. 
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1994). Following the arguments of ZJM03 we used the database of the inclination angles 
compiled by Rankin (1993a,b). This database contains 149 measurements of the magnetic 
inclination angles, whose distribution is presented in the upper panel of Fig. 3. 

3.4. Interpulse emission 

The fraction of interpulses in the observed sample of pulsars provides useful information 
about pulsar geometry (see GH96). However, the statistical studies of this phenomenon are 
difficult since the ratio of amplitudes of the interpulse to the main-pulse is often about 
1% and varies with frequency (Hankins & Fowler 1986). The only representative sample 
of interpulses can be found in Table 6 of the Catalog of 558 pulsars (Taylor, Manchester 
& Lyne 1993). There are 22 pulsars with interpulses (~ 4%) in this sample, with only 
3 certain cases (B0826-34, B0950-I-08, B1929-I-10) corresponding to a single magnetic pole 
(Lyne & Manchester 1988, their Fig. 8). Most hkely PSR B1848+04 with very a broad 
main-pulse also belongs to this category. This corresponds to occurrence rates of about 
3.4% and 0.5% for DP-IP and SP-IP, respectively. These rates were adopted as model values 
by GH96. However, to be consistent with our selection of pulsar periods (section 3.2) one 
should exclude millisecond and other recycled pulsars from both the total pulsar sample and 
from the sample of pulsars with interpulses. When this is done then the corresponding rates 
of occurrence are 2% and 0.8% for DP-IP and SP-IP, respectively. For a comparison, in 
the published sample of 420 pulsars found recently in the Parkes multibeam pulsar survey 
(Manchester et al. 2001; Morris ct al. 2002; Kramer ct al. 2004) one can identify clearly 4 
cases of DP-IP and 1 case of SP-IP, and a few low intensity candidates. This gives lower 
limits of about 1% and 0.24%, respectively. It is difficult to estimate the actual rates of IP 
occurrences until the sensitive search for interpulses in the newly discovered pulsars is made. 



- 12 - 



We assume that the plausible distributions of periods P and inclination angles a in the 
parent pulsar population should be able to explain about 2% of DP-IP and slightly above 
0.5% of SP-IP in the normal pulsar population (excluding millisecond and other recycled 
pulsars) . 

4. Monte Carlo simulations 

We performed the Monte Carlo simulations of the pulse widths Wio, pulsar periods P 
and inclination angles a, in an attempt to reproduce the observed distributions of these 
quantities. We have used the random number generator given in Press et al. (1992). Our 
simulation procedure can be described as a number of subsequent steps: 

(1) Generate the inclination angle a as a random deviate with the parent probability density 

function f{a) corresponding to equations (6), (7), (8) and (9), respectively. 

(2) Generate the observer angle ^ — a + (3 as a. random deviate with the parent probability 

density function /(^) = sin^ (eq. [5]), and calculate the impact angle (3 — ^ — a. 

(3) Generate the pulsar period P as a random deviate with the parent probability density 

function f{P) corresponding to equations (10), (11) and (12), respectively. Record the 
free model parameters of each function used. 

(4) For a given value of P calculate the opening angle pw corresponding to equations (2), 

(3) and (4), respectively. 

(5) Check the detection condition corresponding to equations (13), (14) and (15) for Ntot — 

50000 simulated pulsars for each combination of p — P relations (eqs. [2]- [4]) and 
distribution functions f{P) (eqs. [10]-[12]) and f{a) (eqs. [6]-[9]). Record the number 
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of observed pulsars Nobs- Calculate the beaming fraction / = Nobs/Ntot, as well as the 
rates of detected interpules from two magnetic poles (DP-IP) and a single magnetic 
pole (SP-IP). 

(6) For each set of parameters a, j3 and pio(-P) corresponding to the observable pulsar calcu- 

late the pulse width Wiq according to equation (1), and record the relevant information 
(P^io,P,«). 

(7) Judge the statistical significance using the K-S tests for the simulated and observed 
distributions of pulse width W^q, period P and inclination angle a. 

4.1. Reproduction of previous results of GH96 and ZMJ03 

To make sure that our simulation software works properly, we begun with a reproduc- 
tion of results given in the previous statistical works of GH96 and ZJM03. As expected, 
we managed to reproduce the K-S statistics for Wio, beaming fraction / and the rates of 
observed interpulses given in Table 2 of GH96 (to save space we do not reproduce this table 
here). However, the K-S test for pulsar periods (not included in GH96) resulted in very 
low probabilities (below 10~^), confirming the suggestion of TM98 that the distribution of 
periods in the parent pulsar population is significantly different from that of the observed 
distribution. As we argue later on in this paper, the plausible parent distribution of periods 
can be expressed in the form of a gamma function (eq. [10]) with m — 0.34 ± 0.02 and 
a — 2.52 ± 0.04, in contrast to m = 0.3 and a — 2.5 obtained for the observed sample of 516 
pulsar periods by GH96, or m = 0.28 and a = 2.28 obtained for 1165 periods used in this 
paper (see also ZJM03). 

Twelve entries in Table 1 correspond to cases A1-A4, B1-B4 and C1-C4 from Table 1 in 
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ZJM03 (we excluded cases D1-D4 since for s — 5.8/8.8 = 0.65 they are equivalent to cases 
A1-A4). As one can see, we reproduced quite well the results of K-S tests for P, Wiq and a, 
as well as the values of beaming fraction / (our period sample has few more pulsars below 
0.05 s and above 4.2 s as compared to the one used by ZJM03). We added the rates of the 
interpulse occurrence and concluded that the pure cos a distribution of the parent inclination 
angles (Al, Bl and CI) generates much too few DP-IP (about 0.3% as compared with about 
2% observed). Thus, although the pure cos a distribution gives very good results of K-S 
tests (consistent with the results of TM98 and ZJM03), it should be rejected on the grounds 
of unacceptable interpulse statistics. We found out that the modified cosine distribution of 
ZJM03 (see eq. [9]) not only gives the plausible results of K-S test, but also reproduces the 
rate of occurrence of both DP-IP (~ 2%) and SP-IP (0.5% - 0.7%). 



4.2. Other new results 

We have examined 19200 combinations of distribution functions of pulsar periods (eqs. 
[10], [11] and [12] with 200 combinations of parameters in each case), opening angles (eqs. [2], 
[3], [4a] and [4b]), inclination angles (eqs. [6], [7], [8] and [9]) as well as two options of the 
beam shape: circular (eqs. [14]-[16]) and elliptical (eqs. [17]-[21]). We recorded only those 
cases in which the following conditions were simultaneously satisfied: the probabilities V that 
the observed and simulated distributions of P, Wio and a exceeded 0.1%, as well as the rates 
of occurrence of DP-IP and SP-IP exceeded 2% and 0.5%, respectively. With the adopted 
step of 0.02 in the parameters of gamma, lorentzian and gaussian functions (eq. [10]- [12]), 
this resulted in 15 records. Table 2 presents 5 representative cases with the highest rates of 
SP-IP occurrence. We beheve that case (1) is the most plausible one and Figs. 1-3 present 
a visual comparison of the observed and simulated distribution corresponding to this case. 
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Below we discuss some aspects of our analysis, which are not reflected in Table 2. 

Pulsar period 

Both gamma (eq. [10]) and lorentzian (eq. [11]) functions with parameter values m = 0.34 ± 
0.02, a = 2.52 ± 0.04 and Xq = 0.62 ± 0.02, uq = 0.4 ± 0.02 respectively, are plausible parent 
density distribution functions. However, the gamma function seems much better suited to 
reflect the skew character of the pulsar period distribution. The gaussian function (eq. [12]) 
is rather unlikely. In the best case corresponding to Xq — 0.73 and ctq = 0.36, the significance 
V is only about 0.03% 

Inclination angle 

ZJM03 argued, ignoring the issue of the interpulse emission, that the parent distribution of 
the magnetic inclination angles can be expressed by the cosine function or modified cosine 
function represented by equation (9) in this paper. Including the analysis of interpulse 
statistics we confirm the latter, but we refute the cosine function, since it generates less than 
0.3% of DP-IP (as compared with about 2% observed). This is a strong conclusion since 
better statistics of interpulse occurrences cannot change it in the future. As for other trial 
distribution functions, the sine- like function (eq. [7]) generates too many DP-IP (~ 5%) and 
far too few SP-IP (~ 0.02%), while the interpulse rates generated by the flat distribution 
function (eq. [6]) are roughly comparable with observations, especially when associated with 
the meridionally compressed beam (eq. [17]), but the results of K-S test for the magnetic 
inclination angles are not promising ("P < 10~^). Therefore, the modifled cosine function 
(eq. [9]) is the only plausible density distribution function that satisfies all constraints. We 
have checked whether one can improve the statistical results by changing values of parameters 
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A and B in equation (9). It appeared that without violating the basic condition that V > 
0.1% for all considered quantities (P, Wio and a), one can only increase the interpulse rates 
by a small fraction (e.g. from 0.74% to 0.77% for SP-IP) and thus indeed A ^ 0.6 and 
S ^ 0.15 as suggested by ZJM03. 

Beaming fraction 

For the 15 cases satisfying all constraints adopted in our analysis, the beaming fraction / 
defined as the number of detected pulsars divided by 50000 detection attempts is 0.124 ± 
0.004. This is consistent with the probabilities of observing a normal radio pulsar obtained 
by both TM98 and ZJM03. 

5. Conclusions and Discussion 

In this paper we performed statistical studies using the Monte Carlo simulations of the 
possible parent distributions of pulsar periods P and magnetic inclinations angles a. We gen- 
erated synthetic distributions of the pulse- widths Wio, as well as the interpulse occurrences, 
and confronted them with the observational data by means of the Kolmogorov-Smirnov sig- 
nificance tests. We found out that the observed distributions of pulsar periods, pulse widths 
and inclination angles are relatively easy to reproduce with a variety of trial density distri- 
bution functions. However, when we use the criterion that the K-S significance probabilities 
V for P, 1^10 and a are higher than 0.1% and that the generated interpulse rates agree with 
the observed rates, then we were left with just a few possibilities presented in Table 2. Our 
results can be summarized as follows: 

(1) The distribution of the magnetic inclination angles in the parent pulsar population 
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is a complicated function represented by equation (9), which has a local maximum 
around a = 25° and another weaker one around a = 90°. This function reproduces the 
observed distribution of the observed inclination angles (Fig. 3), as well as generates 
the rates of occurrence of interpulses at the observed levels. 

(2) The parent distribution of periods in the normal pulsar population (excluding millisec- 

ond and other recycled pulsars) can be described by the gamma function (eq. [10]) 
with parameters m — 0.34 ± 0.02 and a — 2.52 ± 0.04. The actual values of these 
parameters differ slightly depending on the p-P relation. 

(3) There is no evidence that the shape of the pulsar beam deviates significantly from a 

circular crossection. It is likely that the internal structure of a typical pulsar beam 
consist of two co-axial cones (eq. [4]) 

(4) The beaming fraction /, that is the fraction of observable pulsars or the probability of 

observing a normal pulsar, is 0.124 ± 0.004. 

In general, our results are consistent with those of ZJM03. However, we added the 
constrains related to the interpulse analysis and found out that this aspect of our analysis is 
the most restrictive one. We were able to reject the pure cosine distribution of the inclination 
angles, but we confirmed their modified cosine function which has a local maximum around 
a — 25° and another weaker one near a — 90°. This means that the evolution of the 
magnetic inclination angles in pulsars cannot be described by any simple law (alignment, 
counter alignment, etc). 

Although we improved the analysis of ZJM03 by adding the interpulse statistics, we are 
still missing an analysis of a possible effect of the intrinsic luminosity of radio pulsars on our 
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results. This problem is, however, very difficult and complicated and we will postpone a full 
treatment to the subsequent paper. The proper approach would be to compare the synthetic 
radio luminosity with the minimum detectable flux achieved in a given pulsar survey, and 
thus it can be applied only to a uniform data sets of pulsars detected in single survey Our 
data do not have such a degree of uniformity However, most surveys were less sensitive to 
long-period pulsars, as it follows from the nature of applied Fourier-transform method. Since 
the interpulse emission (which appears to be the most restrictive constraint in our analysis) 
occur mainly at shorter periods, then a possible underrepresentation of pulsars with longer 
periods should not affect significantly our general results. 

This work is partially supported by the Grant 1 P03D 029 26 of the Polish State Com- 
mittee for Scientific Research. We thank E. Gil and U. Maciejewska for technical help. We 
also thank and anonymous referee for very constructive and helpful criticism. 
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Table 1. Results of Monte Carlo simulations - reproduction of Table 1 in ZJM03 



p-p 


Inclination 


Pulsar 






Kolmogorov- S mirnov 






Beaming 


Interpulse 


relation^ 




angle 


period 






statisticsttt 






fraction 


rates 


Pio(P) 










P 


mo 




a 














m 


a 


V 


V 




V 


V 


V 


f 


DP-IP 


SP-IP 




A. 


1. 


cos a 


0.38 


2.21 


0.049 


0.017 


1.068 


0.236 


0.139 


0.007 


0.119 


0.25% 


0.62% 


A. 


2. 


eq.(9) 


0.38 


2.204 


0.046 


0.029 


0.108 


0.009 


0.102 


0.093 


0.123 


2.24% 


0.36% 


A. 


3. 


2/7r 


0.278 


2.277 


0.151 


10-18 


0.048 


0.666 


0.354 


10-1^ 


0.162 


3.58% 


0.39% 


A. 


4. 


2/7r 


0.411 


2.101 


0.046 


0.028 


0.105 


0.012 


0.360 


10-14 


0.147 


3.37% 


0.41% 


B. 


1. 


cos a 


0.39 


2.339 


0.054 


0.007 


0.121 


0.002 


0.140 


0.006 


0.121 


0.28% 


1.01% 


B. 


2. 


cq.(9) 


0.1 


2.29 


0.051 


0.011 


0.157 


lo-"' 


0.1 


0.110 


0.126 


2.58% 


0.69%i 


B. 




2/7r 


0.27,S 


2.277 


0.181 


0.007 


0.118 


10 1 


0.:-!17 


10 


0.179 


i.-VA 


0.75% 


B. 


4. 


2/7r 


0.4 


2.255 


0.049 


0.017 


0.101 


0.018 


0.354 


10-" 


0.154 


3.72% 


0.49% 


C. 


1. 


COS a 


0.38 


2.27 


0.054 


0.006 


0.086 


0.066 


0.137 


0.008 


0.122 


0.29% 


0.90% 


C. 


2. 


eq.(9) 


0.37 


2.27 


0.057 


0.009 


0.137 


0.019 


0.1 


0.107 


0.128 


2.54% 


0.51% 


C. 


3. 


2/7r 


0.278 


2.277 


0.152 


10-18 


0.07 


0.206 


0.352 


10-14 


0.168 


4.12% 


0.57% 


C. 


4. 


2/7r 


0.378 


2.262 


0.047 


0.023 


0.085 


0.067 


0.355 


10-14 


0.153 


3.57% 


0.43% 





tp - P relation (beam radius): A - pio = 5.°8P-i/3 (eq. [2]), B - pio = 5.°6P-i/2 (eq. [3]), 

C - pio = 4.°9P-i/2 for P < 0.7 s and pio = 6.°3P-i/2 for others (eq. [4b]) 

tt gamma function /(P) = Go(^'/"i)""^e-(^/'") (eq. [10]) 

tttK-S tests: T) — maximum difference between cumulative distribution functions corresponding to the observed and the simulated data samples, 
V - probability that both the observed and the simulated data sets are drawn from the same parent distribution. 



Table 2. Results of Monte Carlo simulations 



p — P Inclination Pulsar Kolmogorov-Smirnov Beaming Interpulse 

No relation angle period statisticst fraction rates 





Pio{P) 








m 
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Ww 




a 
























V 


V 


V 


V 


V 


V 


f 


DP-IP 


SP-IP 




1. 


eq. (4a) 


eq. 


(9) 


m = 


eq. 
0.34 


(10) 
a ■- 


= 2.52 


0.054 


0.006 


0.124 


0.002 


0.093 


0.156 


0.128 


2.07% 


0.74% 


2. 


eq. (3) 


eq. 


(9) 


m = 


eq. 
0.36 


(10) 
a 


= 2.6 


0.055 


0.006 


0.127 


0.001 


0.101 


0.103 


0.121 


2.55% 


0.64% 


3. 


eq. (4b) 


eq. 


(9) 


m = 


eq. 
0.32 


(10) 
a 


= 2.6 


0.054 


0.006 


0.117 


0.004 


0.098 


0.120 


0.125 


2.44% 


0.53% 


4. 


eq. (2) 


eq. 


(9) 


xo = 


eq. 
= 0.6 


(11) 
ao 


= 0.42 


0.053 


0.008 


0.090 


0.049 


0.107 


0.072 


0.123 


2.38% 


0.57% 


5. 


eq. (4b) 


eq. 


(9) 


xo = 


eq. 
0.64 


(11) 
ao 


= 0.38 


0.058 


0.003 


0.115 


0.004 


0.109 


0.061 


0.126 


2.25% 


0.51% 



tsee Table 1 for explanation 
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Table 3: Interpulse emission in normal pulsars (after Table 6 in Taylor et al. 1993) 



No. 


PSR 


P (sec) 


Fractional 


Phase 


How many 








Amplitude 


separation (°) 


poles? 


1. 


B0531+21 


0.033 


0.6 


145 


DP 


2. 


B0823+26 


0.530 


0.005 


180 


DP 


3. 


B0826-34 


1.848 


0.1 


180 


SP 


4. 


B0906-49 


0.106 


0.24 


176 


DP 


5. 


B0950+08 


0.253 


0.012 


210 


SP 


6. 


B1055-52 


0.197 


0.5 


205 


DP 


7. 


B1259-63 


0.047 


0.75 


145 


DP 


8. 


B1702-19 


0.298 


0.15 


180 


DP 


9. 


B1719-37 


0.236 


0.15 


180 


DP 


10. 


B1736-29 


0.322 


0.4 


180 


DP 


11. 


B1822-09 


0.768 


0.05 


185 


DP 


12. 


B1848+04 


0.285 


0.2 


200 


SP 


13. 


B1929+10 


0.226 


0.018 


170 


SP 


14. 


B1944+17 


0.441 


0.005 


175 


DP 



-24- 
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Period (s) 

Fig. 1. — Observed distribution of 1165 pulsar periods with 0.02 s < P < 8.52 s (upper 
panel) and simulated distribution corresponding to case (1) in Table 2 (lower panel). The 
solid lines represent the analytical fit to the observed distribution in the form of a gamma 
function with m — 0.28 and a — 2.28, while the dashed line represents the parent density 
distribution gamma function with m = 0.34 and a = 2.52. 
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Fig. 2. — Observed distribution of 238 pulse widths Wiq taken from GH96 (upper panel) 
and simulated distribution corresponding to case (1) in Table 2 (lower panel). 
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Fig. 3. — Observed distribution of 149 magnetic inclination angles a taken from Rankin 
(1993b) (upper panel) and simulated distribution corresponding to case (1) in Table 2 (lower 
panel) . The solid line represents the parent density distribution function expressed by equa- 
tion (9). 



